Evolution equation of entanglement for general bipartite systems 
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We explore how entanglement of a general bipartite system evolves when one subsystem undergoes 
the action of an arbitrary noisy channel. It is found that the dynamics of entanglement for general 
bipartite systems under the influence of such channel is determined by the channel's action on 
the maximally entangled state, which includes as a special case the results for two-qubit systems 
[Nature Physics 4, 99 (2008)]. In particular, for multi-qubit or qubit-qudit systems, we get a general 
factorization law for evolution equation of entanglement with one qubit being subject to a noisy 
channel. Our results can help the experimental characterization of entanglement dynamics. 
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Introduction — In quantum information theory, entan- 
glement is a vital resource for some practical applications 
such as quantum cryptography, quantum teleportation 
and quantum computation [lj, |2|] . To fulfill such tasks by 
constructing suitable quantum devices, we inevitably en- 
counter some interactions of the multiparticle quantum 
states under consideration with its environment. These 
undesired couplings give rise to decoherence, which de- 
grades the entanglement when the particles propagate or 
the computation evolves. Therefore, it is of great prac- 
tical importance to investigate the dynamics of entan- 
glement for the quantum systems under the influence of 
decoherence. 

Recently much effort has been devoted to understand- 
ing the dynamics of entanglement 0, 0, H, [|| HI • In 
stead of deducing the evolution of entanglement from the 
time evolution of the state, Thomas Konrad et al. 0] pro- 
vided a direct relationship between the initial and final 
entanglement of an arbitrary bipartite state of two qubits 
with one qubit subject to incoherent dynamics, where 
qubit represents the state of 2-dimensional quantum sys- 
tem. It is also discussed in [1] for two-qudit systems with 
either system undergoing an arbitrary physical process, 
where qudit denotes the state of D-dimensional quantum 
system. On the condition that the pure initial state has 
d non-zero Schmidt coefficients, an evolution equality is 
satisfied during an initially finite time interval. 

In fact, for practical applications in quantum infor- 
mation processing, multipartite entanglement are often 
concerned, e.g. cluster states used as a resource for 
one-way quantum computing [9( , multi-photon entangled 
states [lfj etc. In this Letter, we investigate the evolu- 
tion of entanglement for multipartite pure states with 
one part of the system undergoing the action of an ar- 
bitrary noisy channel, which represents the influence of 
an environment, of measurements or of both. Basically 
this is equivalent to studying the evolution of entangle- 
ment for bipartite systems with one subsystem subject 
to some arbitrary noisy channels. Moreover, bipartite 



systems with higher-dimension can improve the perfor- 
mance of various quantum information and computation 
tasks, such as quantum cryptography lf|. Thus it is 
necessary to investigate the dynamics of entanglement 
for an N% x A2 bipartite system under the influence of 
decoherence. In the following, we find that the dynamics 
of entanglement for general bipartite systems under the 
influence of a noisy channel is determined by the chan- 
nel's action on the maximally entangled state, instead of 
exploring the time-dependent action of the channel on all 
initial states. Therefore the robustness of entanglement- 
based quantum information processing protocols is easily 
and fully characterized by a single quantity. As applica- 
tions we discuss two examples in detail: the entanglement 
evolution for a generalized three-qubit W state [12j with 
one qubit undergoing the action of phase noise channel 
and generalized amplitude damping channel respectively, 
and the one for the ground state in a nuclear magnetic 
resonance (NMR) system with one subsystem subject to 
a decay channel. 

Evolution equation of entanglement — We first define 
the entanglement measure for bipartite systems. For a 
pure state \ X ) = E^Ii Ajltf) e cNl ® C ^ 2 in thc 

computational basis |i) and \j) of Hilbert space C Nl and 
C* 2 respectively, we define the concurrence matrix C 
with entries C a p — {x\{L a <8 Lp)\x*), where \x*} is the 
complex conjugate of |x), L a , a = 1,- ■ ■ , N±(N\ — f )/2, 
and Lp, j3 = 1, • • • , A^A^ — l)/2, are the generators of 
SO(Ni) and SO(N2) groups respectively. The Frobenius 
norm of C is just the I-concurrence [13| |. 



C[|X>] = I|C 



F = 



\ 



N 1 (N 1 -l)/2 N 2 (N 2 -l)/2 



E i^i 2 ' « 



p=i 



which reduces to concurrence when restricting to 2 ® 2 
systems 14j . These quantities can be measured for pure 
states |15| . As we konw, the I-concurrence is equal to the 



length of the concurrence vector [H, H3] • 

The Ai ( Ai-f ) 1 2 generators L a of SO(N% ) have the fol- 
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lowing form, L a{kl} =(-l) k+l+1 \k}(l\ + (-l) k + l \l}(k\, k< 
I, where a (fei ) = {ji,h, • ■ ■ ,3n x -2) with 1 < ji < 
j 2 < ■■■ < jNt-2 < Ni, k < I and k,l £ 
{ji, 32, ■ ■ ■ , JNi-z}- Therefore, one has C akl p kn , = 
2{ -l) k + l+k '+ l \A kk ,A w - A kV A lk ,y and C[\ X )} = 
2 \jYln<j Y, k <i \ A ik Aji - A u Aj k \ 2 , which is the just the 
generalized concurrence in [l8| up to a constant factor. 
C[|x)] is zero when |x) is separable, i.e., A^ = biCj for 
some complex numbers b iy Cj. On the other hand, C takes 
its maximum value y/2(N - 1)/N with N = min(iV 1 ,N 2 ), 
when \x) is a maximally entangled state. 

For a bipartite mixed state p — '^2 i Pi\'4 , i)('4 , i\, Pi > 0, 
5^iPi = l, the concurrence is defined by the convex-roof: 



C[p] = mm^piC[\^i 



(2) 



where C[|^)] is the norm of concurrence matrix C[|^)] 
and the minimum is obtained over all possible pure state 
decompositions \ipi) of p. 

Let |x) be the bipartite pure initial state, and let 
the second subsystem undergo the action of a noisy 
channel. We will denote the noisy channel by $ there- 
after. Then the final state of the system takes the form 
p' = (7 (g) $)|x)(x|. To investigate the properties of the 
entanglement of the final state p' , it is convenient to re- 
express the initial state as: |x) = {M x ® I)\4>), where 
M x ^y^E^iEjSiA^Xil is the filtering opera- 
tion [19( | acting on the first subsystem of maximally en- 
tangled pure state, \4>) — J2n=i \ n ) ® \ n )/V^2: an d I 
is the N 2 x N 2 identity matrix. Due to the fact that 
M x and $ act on the first and second subsystem of the 
state \(j>) respectively, the evolution of \x) takes the form: 
pf = {M x <g> I) P$ {Ml <g> /) , where p % = ( I ® $) 1 0) (<f>\ . 

In terms of concurrence matrix, we have the following 
theorem: if p$ is a pure state, the concurrence C[p'] for 
the state p 1 is given by 



CV\= 



\ 



Wi(W 1 -l)/2 

E 

0=1 



N 2 (N 2 -l)/2 

E 

0=1 



(3) 



Hence the concurrence takes the form ([3]). ■ 
Remark: With respect to the relations C[p'} = 
C[\x)]C[p$] in |[ for two-qubit systems (N t = N 2 = 2), 
here we have a similar relation for the corresponding con- 
currence matrices, C[p'] = ^C[|x)]C[p$]. 

For a general channel $, the state p$ is usually a mixed 
one. Assume p$ has an optimal pure state decomposi- 
tion p$ = Ei^l^i)(^l sucn tnat C [Ps\ = EiPi^O^i)]- 

By convexity we have C[(I ® $)|x)(xl] = 
C[J2iPi(M x ® (Mt ® /)] < EiftC[(M x ® 

J)|^i)(^i|(M+ ® /)]. According to the Cauchy in- 
equality, we have | ^^ ( 1 JV2 " 1)/2 C Q7 [|x)]C 7/3 [p$] | 2 < 

E W-D/ 2 | Ca7 [|x)]| 2 E?if 2 - 1)/2 |^[PS]| 2 . In 
terms of Eq. © we get 



C(p') < £>C[(M X ®2)|^)(^|(Mt®2)] 



A/2. 
2 



^c[|x)]E^[l^)] = ^c[l\)/"*^]-(-) 



This inequality can be generalized to the case that the 
initial state pg is a mixed one, 



N 2 



C [(I »$)/*>] < -^C(p )C[p$] 



(6) 



If we consider bipartite states in iVi £g> 2 system, the 
result © can be generalized for arbitrary noisy channels 
$. So we get the following corollary: if the pure initial 
state is a bipartite Ni®2 one, while the second subsystem 
is subject to an arbitrary noisy channel $, we have the 
following evolution equation of concurrence, 



C[p'\ = C[|x)]C[ 



(7) 



This corollary is proved as follows. Without loss of 
generality, we suppose that p$ is a mixed state. By using 
the procedure of the optimal pure state decomposition 
adopted in Ref. [14j, there must exist an optimal pure 
state decomposition for a2®2 mixed state 



= E k '^)^ 



(8) 



where C' a/3 



^E% N2 ~ x)/2 c a ,[\ x )]c y0 [p s ]. <7 Q7 [|x>] 

and C 7/ 3[p$] are the entries of the concurrence matrices 
C(|x)) and C(p$) respectively. 

Let us now prove this theorem. Suppose the pure 
state p$ has the following generic form, p$ = \tp}(tp\, 
= Ei^=i °<iNi)- The nna l state has the form p' = 
where W) = T.^T^l, B M \H\ and B, 

N 2 



/ N 2 E/=i Auaij . Some straightforward algebra yields 
the equation 



C' 



a kil3 k 



ik+l+k'+l' 



2(-l 



„ N 2 (N 2 -l)/2 
1^2 \ ^ 

~2" ^ 

7=1 



{B kk 'B w — B k iiB lk ,y 



such that C[p$] = C[|^i)]; are satisfied. Suppose 
the pure decomposition Eq. (jSJl is not optimal for con- 
currence of p' in terms of Eq. ([2]). Then there must 
exist another decomposition other than Eq. (JSJ, p$ = 
Ei lil^Di^'ili whicb is an optimal pure state decomposi- 
tion of p 1 = (M x <g) I)p$(M x ® I). In terms of Eq. ©, 
we have 



C[p'] = Y,HC[{M X ®IM)] > C[\x)]C[ P% ] 



(9) 



However, in terms of the optimal pure state decomposi- 
tion ([5]) and convexity, we have 



(c akn [\x)}C lPklll [p $ }).(4) C[p']<J2Pi C [(M x ®IM)] = C[\x)]C[p $ ]. (10) 
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It contradicts with Eq. ([9]). Therefore, the optimal pure 
state decomposition (|SJ) is also optimal for concurrence 
® of p'. Therefore, we get Eq. 0. ■ 
The result ([7]) can also be generalized to the case that 
the initial state po is mixed. Let po — ^iP^i/i^l be an 
optimal pure state decomposition, in the sense that the 
average Frobenius norm of the concurrence matrix over 
this pure state decomposition is minimal. According to 
convexity, we have C[(J®$)p ] =C[J2 i Pi(I®$)\'4>i)('>Pi\] < 
EiftC[(J® Using Eq. 0, we have 

C[(I®$)p }<C[p }C{p$]. (11) 

When Ni = 2, the result reduces to the main result 
of Thomas Konrad et al. @]. The results @, @-0, 
and dTTJ) ) show that the dynamics of entanglement for 
bipartite systems under a one-sided noisy channel is de- 
termined by the channel's action on the maximally en- 
tangled state. 

Let us study in which case the evolution equation of 
entanglement holds for 2 ® 2 pure initial state under 
the influence of local two-sided channel $i ® $2. The 
channels $ can usually be expressed as Kraus operators 
[20I ] . First consider the phase noise channel $1 satisfy- 
ing Tr^Si) = Tr(tr„$i) = 0, we have C[($i <g> $ 2 )x] = 
C[($iO$ 2 )|^)((A|]C(|x)) for any initial state |x) = a|00) + 
6|11) in Schmidt expression. If $1 satisfies the condition 
Tr(crJi) = Tr($i) = 0, for \ X ) = a\00) + b\ll), we have 
C[($i®$ 2 )x] = C[($ 1 ®$ 2 )|0)(0|]C(| X )); for | X ) = a|01} + 
6|10), we have C[($i®$ 2 )x] = <7[(J® $ 2 $i)|0)(0|]C(| X )). 

For the case that the initial state is mixed, let us con- 

(a \ 
d* c I • For phase noise channel $ p , 
00/ 

we still have C[(J® $ p )po] = C[p ]C[p$ p ]. 

Application to two realistic systems — Let us consider 
a three-qubit system with the third qubit exposed to 
a phase noise channel $ p . The phase noise chan- 
nel $ p can be expressed as Kraus operators: K\ = 
(0 1 ) > K2 = ( ) ' wn ere the time-dependent Kraus 
matrix elements are v — exp[—Tt] and w = yl — v 2 . We 
study how the residual entanglement [22|,[23| evolves. For 
the initial state |x) J 4sc = a|001)+/3|010)+7|100), general- 
ized W state, it follows that Tr A [(I(g>I(g>$ p )\x)ABc(x\] 

I ® Sj^Tr^ [|x)abc(x|] = {I ® §p)Pbc, since the partial 
trace is a local operation. Therefore, we obtain p B c — 
Tta\x)abc(x\ = |7| 2 |00)(00| + (a|01> +/3|10))«01|a* + 
(10|/3*) and p AC = Tr B | X ) ABC ( X | = |/5| 2 |00)(00| + (a|01)+ 
7|10))((01|a* + (10|7*). Due to the fact that p A B and 
Pac are special cases of po in [2lj, we immediately get 
C[(I <8> % p )pbc] = C[p BC ]C[p $p } = 2|a/3|exp[-n] and 
C[(I ® $ p )p A c] = C[p AC ]C[p $p ] = 2|a 7 |exp[-rt]. If we 
regard AB subsystem as a whole 4 dimensional sys- 
tem, we can calculate C[(I (g> I (g> $ p )/9ab:c] by means 
of Eg 0, C[( 7 «l I ® $ P Ws : c] = C[p AB: c]C[p $p ] - 
2|a|-^|7| 2 + |/?| 2 exp[— Ft]. Therefore the residual entan- 
glement t c{ab) = C[p' AB . c } 2 -C[p' BC ] 2 -C[p' AC ] 2 = 0. 



This indicates that the residual entanglement tc(ab) of 
this initial state keeps invariant when the third subsys- 
tem is exposed to a dephasing channel % p . 

Next we consider the generalized amplitude damping 
channel, $gad, describing the effect of dissipation to 
an environment at finite temperature, which is the case 
relevant to NMR quantum computation. The channel 
usually adopts the form as: K = \/p(l "),Ki = 

Vp(°o Z),k 2 = VT^p~{» l),K 3 = y/T=p(Z 0). 

Without loss of generality, set p = |. In terms of Eq. , 
we obtain C A B-.c(\a\,t) = C[(I <8> I ® $gad)pab-.c] = 
\a\«Jl - |a| 2 (exp[-2rt] + 2exp[-rt] - 1). As shown 
in Fig. 1 (a), over a continuous range of |a| values, 
Cab:C\\ol\,i) actually goes abruptly to zero in a finite 
time and remains zero thereafter. This is the "entangle- 
ment sudden death" (ESD) effect 0,0, El]. The residual 
entanglement is calculated as tc(ab) — C\b c ~ ^ab ~ 
C AC . For simplicity, we focus on \f3\ = \-f\. It is il- 
lustrated in Fig. 1 (b) that the residual entanglement 
terminates completely after a finite interval and remains 
zero thereafter, and it reaches a maximal value at the 
point Tt = 0.0936 and \a\ = 0.4996. 

The evolution of entanglement is basically related to 
the Hamiltonian of a physical system. For example, 
consider the ground state of an NMR system, \ip) = 
■y/T/3\02)- v /T/ 7 3|n) + v / l73~l 2 0), in which the coupling 
Hamiltonian of two spin— 1 nuclei can be expressed as: 




FIG. 1: (a) The concurrence CAB-.c(\a\,i) vs Tt and ampli- 
tude I a I, whree Y is generalized amplitude decay rate. The 
system of AB and of C disentangle completely and abruptly 
in just a finite time for all |q| in the range shown, (b) Depen- 
dence of the residual entanglement Tciab) on Tt and ampli- 
tude |a| for \/3\ — \ The residual entanglement gets maxi- 
mized value at the point Ft = 0.0936 and \a\ = 0.4996. 
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H = JS1S2, J > 0, where Si and §2 are spin oper- 
ators of nucleus 1 and 2 respectively. Supposing that 
is the initial state, under the single-sided relaxation 
operation M = diag(e~ r2t , e~ Tl \ 1), r 2 > T 1 > [2(|, 
the final state becomes p' = (I <g> M)\ip) <g> M^)/p, 
where p = Tr[(7 M)\tp) {4>\{I ® M*)]. In terms of 
Eq. ([3]), we obtain the dependence of entanglement on 
t, c[p'] = y / 4[(e-2r 1 t + i)( e -2r 2f + 1 )_ 1 ]/[ 1 + e -2r lt + 
e -2r 2 tj when i > p / ^ |20). 

Conclusions. — In summary, we have investigated the 
time evolution of entanglement for arbitrary bipartite 
systems, with one part subject to interactions with envi- 
ronments. Explicit expressions are derived for bipartite 
systems and a general factorization law is obtained for 
multi-qubit or qubit-qudit systems with one qubit un- 
dergoing the action of a noisy channel. It allows one to 
know the time evolution of entanglement for an arbitrary 
initial state, if one knows the time evolution of entangle- 
ment for the bipartite maximally entangled state. The 
later only depends on the detailed noisy channel and has 
nothing to do with the initial states. Our results can be 
used to infer the evolution of entanglement under certain 
time-continuous influences of the environment. Due to 
the fact that all the quantities of entanglement measure 
can be evaluated efficiently for pure states, it can help the 
experimental characterization of entanglement dynamics. 
Moreover, the results can be also directly applied to in- 
put/output processes such as gates used in sequential 
quantum computing. As applications we have studied 
the entanglement evolution of the generalized three-qubit 
W state, with one qubit undergoing the action of phase 
noise channel and generalized amplitude damping chan- 
nel respectively. We also obtain the evolution of entan- 
glement for the ground state in an NMR system with one 
subsystem subject to a decay channel. 
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